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Abstract 

We study exceptional modes of both the Wilson and the clover ac- 
tion in order to understand why quenched clover spectroscopy suffers so 
severely from exceptional configurations. We show that, in contrast to the 
case of the Wilson action, a large clover coefficient can make the excep- 
tional modes extremely localized and thus very sensitive to short distance 
fluctuations. We describe a way to optimize the chiral behavior of Wilson- 
type lattice fermion actions by studying their low energy real eigenmodes. 
We find a candidate action, the clover action with fat links with a tuned 
clover term. We present a calculation of spectroscopy and matrix elements 
at Wilson gauge coupling j3 = 5.7. When compared to simulations with 
the standard (nonperturbatively improved) clover action at small lattice 
spacing, the action shows good scaling behavior, with an apparent great 
reduction in the number of exceptional configurations. 



1 Introduction 



There is an increasing accumulation of evidence from lattice simulations of the 
importance of the topological properties of the QCD vacuum (ij |2| . Simulations 
of the pure gauge theory find that the vacuum is filled with instantons of size 
0.3-0.5 fm, with a density of about one per fm -4 , although the precise numbers 
are still controversial ||, An elaborate phenomenology of the properties of 
hadrons is based on instanton dynamics and its connection to chiral symmetry 
breaking Js). The extent to which the complete continuum phenomenology of 
instanton effects in hadronic physics is actually realized by QCD, and can be 
studied using lattice simulations, remains an open question, but its outline is 
present. 

Lattice artifacts associated with instantons also seem to be connected to 
some of the difficulties of numerical simulations with quenched QCD, namely, 
eigenmodes of the Dirac operator which occur away from the critical quark 
mass, and which spoil the calculation of the fcrmion propagator, the so-called 
"exceptional configurations" ||. Precisely how does this occur? And, given that 
instantons are responsible for chiral symmetry breaking, is it possible to opti- 
mize the lattice discretization of a fermion action with respect to its topological 
properties? Our goal is to shrink the range of bare quark mass over which the 
low lying real eigenmodes of the Dirac operator occur, on background gauge 
field configurations which are typical equilibrium configurations at some gauge 
coupling. 

A chiral improvement of the lattice Dirac operator is expected to decrease 
the additive mass renormalization of zero-modes as well as their spread. This 
is indeed the case for classically perfect actions which have been proven to 
have exact zero modes with no shift or spread at all Jjj. On the other hand, 
surprisingly, this does not seem to be true for the simplest improvement on 
the Wilson action, the clover action. Although the clover term reduces the 
additive mass renormalization as compared to the Wilson action, the problem 
with exceptional configurations appears to be more severe for the clover than 
for the Wilson action; at least with the non-perturbatively determined value of 
the clover coefficient C sw On a set of 50 12 3 x24 (3 = 5.9 configurations the 
Fermilab group found 7 exceptional modes both with the Wilson and the clover 
action (C sw = 1.91) but the spread of the clover modes in terms of the pion mass 
was at least twice that of the Wilson modes. This indicates that surprisingly, 
improving the action makes the situation worse for exceptional configurations. 
In the present paper we show why this happens by studying how the real modes 
of the Dirac operator change with the clover coefficient. 

We are then drawn to an alternate action, a modified clover action, where 
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the gauge connections are replaced by APE-blocked || links, and the clover 
coefficient C sw is tuned to optimize chiral properties. Specifically, 

S = y~^(m + 4)i/j(n)^(n) 

n 

- 5 H(^( n )( 1 " 7/x)^NV(n + M) + $(n)(l + l„)V^(n - ^(n - p) 
+ ~Y~ ^( n ) ia »i> p »v{n)i>{n) (1) 

with 

Vi n) (i)= (l-a)^"" 1 ^) 

+yj"- 1) ( 2 ; - ^V^Or - *)K (n_1) (a: -* + £)), (2) 

with V^(x) projected back onto 51/(3) after each step, and (n) = U^(n) 
the original link variable. We take a = 0.45 and N = 10 smearing steps, chosen 
because of our previous work on instantons This choice of parameters is 
not unique and might not even be optimal. P M „ is the usual clover set of links, 
but built of the V^'s. At Wilson gauge coupling 5.7 — 5.8, the best choice is 
C sw = 1.2, and it decreases to the tree-level C sw — 1 value at larger (3. We will 
also quote simulation data in terms of the hopping parameter, k = 1/(8 + 2m). 

The optimized action has a spread of low lying real eigenmodes with respect 
to the bare quark mass which is less than a third of the spread for the usual 
Wilson action. In terms of the square of the pion mass the spread is about three 
times smaller for the optimized action than for the standard Wilson action. The 
action has other good features, as well: the renormalization factors connecting 
lattice quantities to their continuum values appear to be very close to unity. 
The action also appears to require only about half the number of sparse matrix 
inversion steps as the usual clover action (at equivalent values of the physical 
parameters) . 

The action we propose is completely unimproved in its kinetic properties, 
so its dispersion relation and heavy quark mass artifacts are identical to those 
of the Wilson action. It would be very easy to improve it by beginning with a 
more complicated free fermion discretization. 

Fermion actions with fat links have received considerable attention in the 
past year. The first use of them we know of (although with a different motiva- 
tion) was in the simulations of QCD on cooled gauge fields by the MIT group 
[ p"o| . The MILC collaboration (T^|, Orginos and Toussaint Jl2| and Lagae and 
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D. K. Sinclair [|T3| have shown that modest fattening considerably improves fla- 
vor symmetry restoration in simulations with staggered fermions. We have used 
calculations of staggered spectroscopy on highly smoothed gauge configurations 
to compare chiral symmetry breaking in SU (2) gauge theory and in instanton 
backgrounds |L4|]. Finally, all fixed point actions and approximate fixed point 
actions we know of |l^, [16| for fermions seem to incorporate fat links. Fixed 
point fermions realize the index theorem and retain chiral symmetry at nonzero 
lattice spacing |?j , and so one way of viewing a fat link action is as an approxi- 
mate FP action, which includes its chiral properties but not its kinetic ones. 

An apparent drawback of fat link actions is the lack of a transfer matrix 
between consecutive time slices. N APE steps can mix links up to ±N timeslices 
away and a strict transfer matrix cannot be defined on time slices closer than 
2N lattice spacings (20 in our case). In practice we are concerned only with 
the exponential decay of correlation functions and we found asymptotic decay 
after 2-5 time slices in our simulation. This is expected if we realize that APE 
smearing is basically a random walk whose range can be estimated as \/N a « 1.4 
in our case. 

The outline of the paper is as follows: In Section 2 we review the continuum 
index theorem and describe how lattice fermions fail to reproduce it. We then 
describe how we find the real eigenmodes of the Wilson-Dirac operator. All 
this is basically a review, and experts may skip it. In Section 3 we describe the 
connection between small topological objects and exceptional configurations for 
the thin-link clover action. In Section 4 we return to the fat link clover action 
and describe the tests we performed to tune the action. Section 5 describes a 
calculation of spectroscopy and matrix elements at aT c = 1/4 (Wilson gauge 
coupling = 5.7) using the new action. 

2 Measuring the Chiral Properties of Wilson- 
like Fermions 

The reader might recall that in the continuum the local topological density 

q( x ) = ttt; ^t-iivpoF^uFpc (3) 
oZir 

is related to the divergence of the flavor singlet (with n f flavors) axial- vector 
current 

= 2mTpi^ijj + 2injq. (4) 
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Integrating over all x, this relation implies a connection between the topological 
charge Q — J d 4 xq(x) and a mode sum, 

Q = mTr^jsip = mTr^s -J- = mV -^Ts/s (5) 

7 • L> + m iA s + m 

where f s are the eigenfunctions of the (antihermitian) Dirac operator 7 • D, 
7 ■ Df s — iX s fs, with A s real. The property {75,7 ■ D} = leads to the 
condition that fj^fs = if A s 7^ 0, and if A s = we may choose f s to have a 
definite chirality, 75/3 = ±/ s . Thus it follows that 

Q = = /sTs/s = »+ - n - ( 6 ) 

s,\ g =0 

where n + and n_ are the number of zero eigenmodes with positive and negative 
chirality. This is the index theorem. 

On the lattice, essentially every statement in the preceding paragraph is 
contaminated by lattice artifacts p7| . Here we focus on Wilson-like fermions, 
where chiral symmetry is broken by the addition of terms proportional to the 
Dirac scalar and/or tensor operators. The Wilson or clover fermion action 
analog of 7 • D, D w , is neither Hermitian nor antihermitian and its eigenvalues 
are generally complex. D w can also have real eigenvalues. These real eigenvalues 
usually do not occur at zero. Their locations spread around — m c where m c is 
the critical bare mass, at which the pion mass vanishes. For Wilson-like fermions 
m c is usually negative and its magnitude decreases with increasing clover term. 

On smooth, isolated instanton background configurations, the location of the 
real eigenmode varies with the size p of the instanton. For large instantons, the 
low lying real mode occurs close to zero. Accompanying this near-zero mode, 
there is a set of modes which do not cluster around zero, but around ~ O(l) 
(in lattice units). These are "doubler modes." As one decreases the instanton 
size, the eigenmode near zero shifts in the positive direction, approaching the 
doubler modes. This p-dependent mass shift and the location of the doubler 
modes all depend on the particular choice of lattice fermion action Jl8|, |l9) . The 
shift of the low lying real eigenmode for a given instanton size for the Wilson 
action is larger than for the clover action, which has better chiral properties. 
Nevertheless, for both actions, as the instanton size decreases, sooner or later 
the lowest eigenvalue starts to increase, approaches the doublers and eventually 
annihilates with one of the doublers - the fermion does not see the instanton 
any longer. 

Real eigenmodes of the Dirac operator also occur on equilibrium background 
configurations. Typically, they occur around some average value, an overall 
(3— dependent mass shift. The real eigenmodes also sense the topological content 
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of the background configuration and their positions can sometimes be correlated 
with the sizes of the background instantons. This is illustrated in Fig. [I] for 
the Wilson action, where we show the instanton size versus the corresponding 
fermionic eigenvalue both for smooth instantons and instantons identified on 
real Monte Carlo generated configurations at j3 = 6.0 on 12 4 lattices. The 
instanton sizes were measured using the method of Rcf. || and the real modes 
were extracted as described below. To find the corresponding instanton for 
each real mode, we compared the quark density of the given mode to the profile 
of the charge density. Usually the quark density had a well identifiable lump 
sitting exactly on an (anti)instanton. These lattices are fairly small in physical 
units with volumes V ~ 2 fm and contain on average about two topological 
objects. In some cases, however, these modes were peaked at not only one 
instanton, but more, and sometimes the quark density was so spread out that 
it was hard to identify a well defined peak. Fig. [j] contains results only from 
uniquely identifiable modes. 

In addition to the low lying modes there are also doublers. Again there are 
two issues: First, do the low lying modes separate from the doubler modes, and 
second, what is the spread of the low lying modes. 

The first question is important because, if the low lying modes do not sepa- 
rate from the doubler modes, then the physics of chiral symmetry breaking with 
a lattice cutoff is different than in the continuum. The lattice theory is no longer 
an approximation to the continuum, and one cannot speak of chiral symmetry 
breaking as being induced by instantons. Chiral symmetry is certainly broken 
in QCD for any value of the cutoff, including the strong coupling limit pl| , but 
the mechanism does not involve instantons. 

Note that if the low lying modes do not separate from the doubler modes, it 
does not make sense to talk about the spread of the low lying modes either. 

The determination of the real eigenvalues A of D w is equivalent to computing 
the bare mass values m\ = —A, where the quark propagator {D w +m) becomes 
singular. In the remainder of the present paper instead of the eigenvalue, we 
shall always quote the corresponding bare mass value, where the propagator is 
singular. On any particular background gauge configuration this can be com- 
puted very simply. We approximate 



with a noisy estimator: cast a random vector rji on each site i, solve (D w +m)x = 
rj for x and measure 



P(m) = (V'75'0) = Tr(Au + m) Ss 



(7) 




(8) 



G 




Figure 1: The instanton size (in units of the lattice spacing) vs. the eigenvalue 
of the corresponding (Wilson action) fcrmionic real mode on smooth instantons 
(crosses) and on real Monte Carlo generated configurations at (3 = 6.0 (squares). 
The horizontal line indicates — m c for the j3 — 6.0 quenched ensemble. 



as a function of to. If an eigenmode of D w is located at A, its presence is sig- 
naled by the appearance of a pole in P(m) at m = —A. This is a variant on 
the standard method of measuring (■tpt/j) in a lattice simulation. It has also 
seen considerable use by the Fermilab group |j) in their studies of exceptional 
configurations (they use a flat source r)i — constant, not a noisy source). The 
eigenfunction itself of a particular mode can be found by performing the inver- 
sion for a test mass very close to the given pole. Then (D + m) _1 projects out 
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the corresponding eigenmode from the source. In this case XX 1S (proportional 
to) the probability density for that mode. For computing the probability density 
we used a flat source rather than a random one. 



3 Trouble with the Thin Link Clover Action 

As we have seen in Fig. |], on fine enough lattices with the Wilson action there 
is a strong correlation between the instanton size and the eigenvalue of the cor- 
responding fermionic real mode; (near) exceptional modes correspond to large 
instantons. We shall now explore what happens in the presence of the clover 
term. We start by looking at smooth single instanton configurations. Our 
naive expectation is that since the clover term is designed to improve the chi- 
ral properties of the action, its inclusion should make the fermion modes more 
"continuum-like" . In other words, it should move the modes corresponding to 
different instanton sizes closer to one another and also shift them all towards 
zero. This is indeed what happens. In Fig. || we plot how the fermion eigen- 
values associated with instantons of various sizes change as a function of the 
clover coefficient. For better legibility in the figure we include only the range of 
C sw > 1.0. For the Wilson action (C sw = 0) the eigenvalues corresponding to 
instantons of size 1.2 < p < 2.5 are spread between 0.19 and 0.86. 

As Csw increases from zero to the tree level value, C sw — 1, all the modes 
move closer to zero and their spread decreases. If C sw is further increased, 
the modes start to pass through zero and also the trajectories corresponding to 
instantons of different sizes start to cross one another. By C sw — 1.4 the order 
of all the eigenvalues for instanton sizes > 1.5a has been reversed. As C sw is 
increased, modes of smaller and smaller instantons cross over to the other side 
of the distribution. 

What happens on equilibrium Monte Carlo configurations? We first per- 
formed a detailed study of the Fermilab quenched j3 = 5.9 exceptional configu- 
rations reported in Ref. ^ . We determined the wave function of both the seven 
Wilson and the seven clover exceptional modes together with the instanton con- 
tent of the corresponding background gauge configurations. Whereas in most of 
the cases the Wilson modes were sitting on a single large instanton, the clover 
modes were always spread over several topological objects. 

This suggests that the mechanism responsible for the exceptional configura- 
tions might be different for the Wilson and the clover action. If indeed there is 
such a difference, we expect it to be more pronounced on coarser lattices, where 
exceptional configurations are more of a problem and also the clover coefficient 
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Figure 2: The real fcrmionic eigenvalue of the clover action versus the clover 
coefficient on instantons of sizes p/a = 1.2,1.4,1.6,2.0,2.5. Bigger symbols 
correspond to larger instantons. 



is usually set to larger values. To explore this further, we studied an ensemble 
of P = 5.7 6 3 x 16 lattices. 

We followed how the wave function of a few typical very exceptional modes 
changed with the clover coefficient starting from the Wilson action (C sw — 0.0) 
up to the non-perturbatively determined value (C sw — 2.25). Ideally, we would 
have liked to compare the peaks of the quark density in the wave function with 
the instanton locations and sizes. Unfortunately, at this rather large value of 
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eigenvalue 


location 


width (a) 


tp^ip max/min 


0.999 


2 10 6 


2.8 


1.8xl0 2 


1.041 


4 2 13 


3.3 


1.4xl0 2 


1.241 


5 5 


2.3 


5.3xl0 3 


1.530 


2 1 2 14 


2.1 


2.1xl0 5 


1.610 


5 1 


1.4 


5.8xl0 5 



Table 1: The physical real eigenmodes of the Wilson Dirac operator (C sw = 0.0) 
on a typical (3 — 5.7 6 3 x 16 lattice. We also show the location where the quark 
density is peaked, the width of the peak (in lattice units) and the ratio of the 
maximum and minimum quark density in the given mode over the whole lattice. 

the lattice spacing the instantons are quite small and the fluctuations large, so 
the size of the instantons is not well defined. We can still ask however how the 
fermions with the Wilson action "feel" whatever is left of the topology on these 
coarse lattices. We have therefore computed the wave function of the real modes 
in the physical branch with the Wilson action. For each wave function we also 
determined the position of the largest peak and its width. Usually the magni- 
tude of the second largest peak was much smaller. The precise definition of the 
width is irrelevant here. Since we use the width only to qualitatively compare 
the localization of different wave functions, any sensible definition would give es- 
sentially the same results. We summarize the data for one typical configuration 
in Table 0. 

The critical mass is m c = —1.04, so the first mode is exceptional. The last 
mode in the Table is approximately halfway towards the doublers. There is a 
clear tendency that the quark density in modes close to —m c has an extended 
broad peak and for modes towards the doublers the peaks get narrower and 
sharper. This is in complete agreement with |^(|. Let us now follow the con- 
tinuous evolution of the Wilson mode at 0.999 if the clover term is gradually 
turned on. Up to C sw — 1.0 the main peak of the quark density remains at 
the same location. It only slowly gets narrower: at C sw — 1.0 its width is 
2.3. Increasing C sw further, another peak appears and its relative significance 
increases with the clover coefficient. By C sw = 2.25 the wave function is com- 
pletely concentrated on a very sharp peak (of width w 1) located at (0 5 1). 
This is the location of the A = 1.610 Wilson-mode which is halfway between the 
doublers and —m c . The wave function of this C sw = 2.25 exceptional mode is 
very similar to that of a Wilson mode lying halfway between the physical and 
the doubler branch. 

We can now easily put the picture together and describe qualitatively how 
the exceptional modes change with the clover coefficient. The first Wilson mode 
can be roughly thought of as being concentrated on a single large instanton (or 
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instantons with roughly equal size). As the clover term is turned on, the zero 
modes corresponding to different instanton sizes get closer and the quark wave 
function spreads over several topological objects, as we also saw on the Fermilab 
configurations. If C sw is further increased, the zero modes separate again, but 
this time the well localized modes (corresponding to smaller instantons) have 
smaller eigenvalues. The relative importance of the broader peaks decreases 
and the mode can become entirely concentrated on a very sharp peak. We 
want to emphasize however that on a given ensemble of gauge configurations at 
a given value of C sw different (near) exceptional modes can look qualitatively 
very different. Some of the modes are entirely concentrated on a sharp peak 
— typically these are the most exceptional ones — some have several peaks of 
various widths. The exact picture presumably depends on the arrangement of 
topological objects and fluctuations on the given configuration. 

The sharply peaked modes look very much like (near) doubler modes of the 
Wilson action and the corresponding eigenvalues and the way they change with 
the clover term can be very sensitive to the fluctuations on the shortest distance 
scale. We indeed encountered several exceptional eigenvalues at C sw = 2.25 
(the non-perturbatively determined value) that moved anomalously quickly with 
the clover coefficient. This is the reason why the clover coefficient cannot be 
optimized to minimize the spread of the real eigenvalues. 



4 Testing and Tuning Actions 

Our goal is still to tune a clover action to have good chiral properties, i.e. to 
minimize the spread of the low lying modes. We see that if the clover coefficient 
is made too large, that program will fail. We also know from previous work [jl4j 
that an action with fat links is insensitive to short distance fluctuations, but 
still knows about instantons and the additional long distance behavior of the 
gauge field responsible for confinement. So we will begin with a fat link action 
of Eqn. [j]. Throughout this paper we create the fat link by 10 APE smearing 
steps with smearing coefficient a = 0.45. 

Fig. H shows the distribution of the real eigenmodes at j3 — 5.8 in terms 
of the bare quark mass, where the propagator is singular, for clover coefficients 
C sw — 0.0 (Wilson action), C sw = 1.0, C sw = 1.2 and C sw = 1.4. The poles 
were located using the pseudoscalar density function of Eqn. || on 40 8 4 config- 
urations. The horizontal ranges of the histograms are equal to the range the 
eigenmode search had been performed. To the right of the distributions is the 
confining phase. If one measures the pion mass in the positive mass region and 
extrapolates to zero with the bare mass mo, one finds that the bare mass at 
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which 771 tt = lies at a value m c located within the range of the low-lying real 
eigenmodes. If a particular configuration has an eigenmode with m\ > m c , the 
quark propagator will be singular here. To the left of the distributions, in the 
negative quark mass region, are the doublers. In Fig. |^ there is no indication of 
them, since the doublers did not show up within the range we scanned for eigen- 
modes. It appears that the doublers and low lying modes are well separated at 
(3 = 5.8. 

Using several fermion actions — all different from the one in the present 
paper — the authors of Ref. [|o| found that real modes cover the whole investi- 
gated negative quark mass range, for all values of gauge couplings they studied, 
as long as the volume is large enough. Our conclusions are not in contradiction 
with this result. They are much more modest-that the action can be tuned 
to make the distribution of real eigenmodes sharply peaked near a small mean 
value, and that the doublers are well separated. No intrinsic property of our 
action protects chiral symmetry at nonzero lattice spacing and so in principle 
it should have real eigenmodes anywhere as well. But for practical calculations 
there is a big difference between a distribution with narrow peaks or a uniform 
spread of eigenmodes, and we are trying to construct the first kind of action. 

To quantify the spread of the low energy real modes, in Fig. |^a we plot the 
average real eigenmode location as the function of the clover coefficient. The 
error bars here are not errors, they are the spread of the modes in m. Fig. 
shows only the spread as the function of the clover coefficient. m c , where the 
pion becomes massless, lies somewhere in the middle of the mode distribution, 
closer to its top (large mass) end. The average eigenmode locations and the 
upper end of the error bars in Fig. bracket m c . Smearing the link removes 
most of the additive mass renormalization even for the Wilson action, changing 
m c ~ —0.95 for the thin link action to m c ~ —0.22 for our case. Adding 
a clover term to the action further reduces the additive mass renormalization, 
and even larger clover terms induce a positive mass renormalization. At f3 — 5.8 
the additive mass renormalization is minimal for c » 1.2. This is also the value 
where the spread of the eigenmodes is minimal. To conclude, the results indicate 
that at (3 = 5.8 (lattice spacing a ~ 0.15 fm) the low lying and doubler modes are 
well separated and explicit chiral symmetry breaking is minimized with clover 
coefficient C sw = 1.2 with our fat link action. 

The situation is less convincing at (3 = 5.7 (lattice spacing a ~ 0.20 fm). Fig. 
||, again based on 40 8 4 configurations, shows the distribution of the eigenmodes 
at (3 = 5.7, clover coefficient C sw — 1.0, 1.2 and 1.4. Here the low lying modes 
are not as well separated as for (3 — 5.8. For C sw = 1.0 the distribution has 
a large tail extending towards negative quark masses. The situation is a bit 
better for C sw — 1.2. Since the low lying modes are not well separated from 
the doublers, it makes no sense to calculate the spread of the distribution. The 
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continuum description of chiral symmetry breaking is about to break down at 
j3 = 5.7. If any of the distributions of Fig. || describes continuum physics, it is 
C sw = 1.2 where the overall mass renormalization is close to zero and the low 
lying modes are best separated form the doublers. 

Finally, at f3 = 5.55 (lattice spacing a ~ 0.24 fm) Fig. |] shows that the 
distribution of eigenmodes is broad and the low lying modes and doublers are 
completely mixed. 

These figures show that it is not possible to make the lattice spacing greater 
than about 0.2 fm, and still retain the continuum-like description of chiral sym- 
metry breaking using a clover-like action. Other, better tuned actions might 
perform better at large lattice spacing. Our exploratory studies with a hyper- 
cubic fixed point jl6| action showed that the FP action is not better in terms 
of chiral symmetry breaking than the fat link clover action with C sw = 1.0. Its 
dispersion relation was improved at lattice spacing a = 0.36 fm, and it (and 
the clover action) showed only a small amount of scale violation in hyperfine 
splittings at large lattice spacing. This just shows that scaling or near scaling of 
a few quantities does not guarantee that the physical mechanism responsible for 
chiral symmetry breaking is the same as in the continuum. This should serve 
as a warning sign for any calculations at large lattice spacing with actions of 
untested chiral properties. 

In principle, one should optimize the spread not in terms of the bare quark 
mass, but in terms of some corresponding physical observable, like the pion 
mass. However, at large fattening, the relation between the bare quark mass 
and the pion mass shows little variation with C sw . 

The usual nonperturbative tuning of the clover term p2f attempts to opti- 
mize current algebra via PCAC relations. It is not clear to us how this is related 
to what we do, since it has been done so far only for thin link actions. We know 
that the optimal C sw by our criterion for a thin link action is larger than the 
usual nonperturbative tuning, for example 2.5 vs 2.0 at /3 = 5.8. Our criterion 
is based on a direct attack on low energy real eigenmodes. To the extent that 
the dynamics of QCD for light quark masses is dominated by this physics, we 
feel our optimization criterion is well founded. 

We do not know if there is a connection between this tuning method and 
the Ginsparg- Wilson J2^| realization of chiral symmetry on the lattice. The free 
field limit of these fat link actions is identical to the standard Wilson action, 
and so while we arc improving the chiral properties of the interacting theory, 
we did not improve the chiral properties of the free theory. 
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Figure 3: The distribution of the bare masses m, where the quark propagator 
is singular at (3 = 5.8 of fat link Wilson (a); fat link clover fermions C sw = 1.0 
(b); C sw = 1.2 (c); C sw = 1.4 (d). 



14 



0.1 



0.07 



<> - 




0.0 



I 



<> 



0.1 



B 



0.2 



o 



0.3 



0.4 



0.02 



0.0 



0.5 



1.0 



1.5 



0.0 



0.5 



1.0 



1.5 



C 



c 



Figure 4: (a) The average real eigenmode location as the function of the clover 
coefficient at (3 — 5.8. The error bars show the spread of the modes, (b) The 
spread of the modes from (a) as the function of the clover coefficient. 

5 Numerical Tests 

We have done a calculation of spectroscopy and matrix elements at lattice spac- 
ing aT c — 1/4 ((3 ~ 5.7 with the Wilson gauge action) using the fat link clover 
action with C sw = 1.2, on a set of 80 8 3 x 24 Wilson f3 = 5.7 configurations. 
This is on the edge of the allowed range of lattice spacings according to our 
criterion of the last section. By itself, one lattice spacing is not a scaling test, 
but we can combine our results with those from other actions to compare the 
new action to them. 

We expect that minimizing the spread of real eigenvalues will also improve 
the situation with exceptional configurations. This is indeed what happens. As 
a preliminary calculation, we determined all the exceptional modes on a set of 
40 6 3 x 16 Wilson (3 = 5.7 configurations both for our optimized fat link action 
and the thin link clover action with C sw — 2.25. The histograms of Fig. fj] show 
where the exceptional modes occur in terms of the lattice pion masses for both 
actions. For the fat link action all the modes occurred at pion masses m, < 0.3 
which corresponds to m v /m p < 0.38, whereas the thin link clover action had 
five exceptional modes above = 0.3 at = 0.42, 0.53, 0.39, 0.30, 0.62 (in 
lattice units). 

On the 8 3 x 24 spectroscopy data set we only identified the eigenmodes 
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Figure 5: The distribution of the bare masses m, where the quark propagator 
is singular at (3 = 5.7 of fat link clover fermions C sw — 1.0 (a); C sw = 1.2 (b); 
C sw = 1.4 (c). 
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Figure 6: The distribution of the bare masses to, where the quark propagator 
is singular at j3 = 5.55 of fat link clover fermions C sw = 1.0 (a); C sw = 1.2 (b). 



occurring above m v = 0.3, as going any closer to the critical mass would have 
been prohibitively costly. With the thin link C sw — 2.25 action we found 11 
poles above = 0.3, the most outlying being at = 0.7. Using the fat link 
action no pole was found above = 0.3. 

We also noticed that the new action seems to be more convergent than the 
thin link clover action: the same biconjugate gradient code needs about half as 
many steps to converge to the same residue as the usual thin link clover action, 
for the same n/p mass ratio. 



5.1 Spectroscopy 

The spectroscopy measurement is entirely straightforward. We gauge fixed to 
Coulomb gauge and used a Gaussian independent particle source wave function 
ijj(r) = exp(— (r/ro) 2 ) with ro = 2. We used pointlike sinks projected onto low 
momentum states. We used naive currents {^^ip, etc.) for interpolating fields. 
The spectra appeared to be asymptotic (as shown by good (correlated) fits to a 
single exponential) beginning at t ~ 3 — 5 and the best fits were selected using 
the HEMCGC criterion ||. 

Our fiducials for comparison are Wilson action and clover action quenched 
spectroscopy. We have tried to restrict the data we used for comparison to 
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Figure 7: The distribution of exceptional fcrmionic modes for the thin link clover 
action with the non-perturbatively determined value of C sw = 2.25 (dashed line) 
and the optimized fat link action (full line) . The eigenvalues are quoted in terms 
of the pion mass (in lattice units) where the corresponding eigenmode makes 
the quark propagator singular. 



lattices with the proper physical volume. 



We can roughly estimate the critical bare quark mass (at which the pion is 
massless) by linearly extrapolating m \ to zero in m n . We also estimated the 
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Figure 8: Bare squared pion mass (crosses) and quark mass from Eq. 12 (di- 
amonds) vs hopping parameter for the C aw =1.2 action. The octagon shows 



critical bare mass using the PCAC relation 

V p ■ (Vn'sV'WVn'sTMV'^)) = 2m q (ipj5ip(0)ipj5ip(x)). (9) 

which, going to the lattice and following |25| , is done by fitting the pseudoscalar 
source-pseudoscalar sink to 

P(t) = Z(exp(-m 7r <) + exp(-m w (N t - t))) (10) 

and the pseudoscalar source-axial sink to 

= _P^9 Z(exp( _ (m ^ _ exp („ m ^ (i v t _ ( U ) 
Z A 

to extract m q . Fig. || shows the squared pion mass vs bare quark mass. The 
quark masses from the local pseudoscalar and axial currents are also shown. 

Fits of the squared pion mass (am,) 2 = B(1/k — 1/k c ) give B — 1.55(2) 
and k c = 0.12515(7). The quark mass is fit to am q = A(1/k — 1/k c ) and 
we find A = 0.447(7), k c = 0.12515(6). In free field theory, we would expect 
A = 1/2, k c = 1/8, and this is our first hint that perturbative corrections to 
bare quantities are small. 

As a scaling test we compare m p /T c and mjv/T c vs. m n /T c for the C sw = 1.2 
action in Figs. || and |l^. 
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Figure 9: Octagons show m p /T c vs. m^/T c for the C sw = 1.2 action at aT c 



1/4. Also shown are Wilson action data, with crosses for aT c 
aT c = 1/8, and fancy crosses for aT c = 1/12. 
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Figure 10: «!» /T c vs. m^/Tc for the C st( , = 1.2 action and the Wilson action, 
labeled as in Fig. 6. 
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K 


PS 


V 


N 


A 


0.118 
0.120 
0.122 
0.123 


0.868( 3) 
0.730( 4) 
0.567( 5) 
0.467( 6) 


1.053( 6) 
0.965( 8) 
0.880(12) 
0.837(17) 


1.589(13) 
1.415(19) 
1.219(28) 
1.093(36) 


1.725(19) 
1.593(22) 
1.501(24) 
1.432(27) 



Table 2: Table of best-fit masses, C S w — 1-2, p — 5.7. 



We compare scaling violations in hyperfine splittings by interpolating our 
data to fixed ir/p mass ratios and plotting the N/p mass ratio vs. m p a. We do 
this at three n/p mass ratios, 0.80, 0.70 and 0.60, in Fig. [H]. In these figures the 
diamonds are Wilson action data in lattices of fixed physical size, 8 3 at j3 = 5.7 
[H), 16 3 at p = 6.0 § 24 3 at (3 = 6.3 §7|, and the crosses are data in various 
larger lattices: 16 3 and 24 3 at = 5.7 and 32 3 at /3 = 6.17 ||, 24 3 at (3 = 6.0 
[|J. When they are present the data points from larger lattices illustrate the 
danger of performing scaling tests with data from different volumes. The bursts 
are from the nonperturbatively improved clover action of Refs. and p9|| . 
The squares show the C sw = 1.2 action. 

We give a list of masses from the C sw = 1.2 action in Table |[ 

We conclude that the fat link clover action has at least as good scaling 
behavior as the usual nonperturbatively tuned clover action. 



5.2 Dispersion Relations 



To view the dispersion relation, we first plot E(j>), the energy of the state 
produced with spatial momentum p, as a function of \p\. The result for the 
C sw = 1.2 action at k = 0.118 is compared to the free dispersion relation at 
aT c = 1/4 in Fig. [l^. It looks very similar to results from the Wilson action at 
the same lattice spacing-an entirely unsurprising result. 

This behavior is quantified by measuring the squared speed of light, c 2 = 
(E(p) 2 —m 2 )/p 2 , iorp = (1, 0, 0). We do this by performing a correlated fit to the 
two propagators. The result is presented in Fig. |l3| and shows the worsening of 
the dispersion relation at larger quark mass, characteristic Wilson/clover kinetic 
behavior. 
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Figure 11: A scaling test for the C sw = 1.2 action (square) vs. Wilson actions 
on lattices of fixed physical size (diamonds) and larger volumes (crosses), and 
the nonperturbatively improved clover action (bursts). Data are interpolated 
to ir/p = 0.80 (a), 0.70 (b), and 0.60 (c). 
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Figure 12: Dispersion relation for hadrons at aT c = 1/4 (a ~ 0.18 fm) from 
the C sw = 1.2 action. The curves are the continuum dispersion relation for the 
appropriate (measured) hadron mass. 

5.3 Renormalization factors 

We now turn to a (rather naive) set of measurements of simple matrix elements. 

Formally, the fat link action can be regarded as "just another" 0(a 2 ) im- 
proved action, if one assumes that the coefficient of the clover term is actually 
C sw = 1 + 0(g 2 ). In perturbation theory, which involves the vector potential 

rather than the link, the action shows the usual cancellation between the 
clover term and the scalar -ipipA term. The spectrum is 0(a 2 ); matrix elements 
using "rotated fields" ip — ► (1 — (7 • D)/2)ip are 0{a 2 ) improved. Of course, the 
lattice-to-continuum renormalization Z-factors are different than the usual thin 
link clover Z-factors. 

In principle, one could calculate all Z-factors using lattice perturbation the- 
ory. The Feynman rules differ from the rules for the usual clover action, only 
in that the vertices are multiplied by form factors, which are functions of the 
gluonic momentum. We have not tried to do this yet. Since we are working at 
fairly strong coupling, we felt that it would be easier to compute the Z-factors 
nonperturbatively, beginning with the vector current and then computing the 
axial renormalization factor using Ward identities. Since this is just the first 
calculation using this action, we will restrict ourselves to the naive (local) cur- 
rents. 
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Figure 13: Squared speed of light vs. hadron mass in units of T c , for (a) 
pseudoscalars, (b) vectors) and (c) protons, from the C sw = 1.2 action at (3 = 
5.7. 
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We begin with the vector current. The conserved (Noether) current for the 
fat link action is the same as the Wilson current, just with fat links instead of 
thin links, 

Jt onS i n ) = i(^N(7A l -l)^W(^ + A)-^(^ + A)(7M + l)^ t N^(«))- (12) 

This current is conserved but not improved. We choose to define the vector 
current Z-factor from the ratio of forward matrix elements 

v _ (Tr(T)Tr(O)) 



v (ir(T)J (t)ir(0)) 

for < t < T. The local current is defined as Jj? c (n) = ip(n)-fni[>(n). We 
measured the Z-factors using 20 configurations on a periodic 8 3 x 24 lattice 
with T = 10 and averaged 9 > t > 1. The results of a jacknife analysis are 
shown in Table ||. 



The fact that Z cons ^ 1 is a finite size effect. In the numerator of Eqn. 13 
both quarks can propagate in both time directions around the torus, but the 
forwards and backwards paths contribute differently to the denominator. We 
correct this phenomenologically by computing the ratio z loc /Z cons and display 
it in the table. 



Lattice perturbation theory which correctly takes into account the residue 
of the pole of the massive quark [^0) predicts that this ratio is equal to 

Zv = ^-^Z V (14) 

where the (1 — 6k) factor is the residue (recall k c = .125) and Zy = l + ag 2 + . . .. 
As Table y shows, the correction to the tree level formula, Z, which is what is 
usually quoted as the Z-factor, differs from unity by less than than two per cent 
at k = 0.123. 



We also attempted to measure the axial vector current renormalization using 
the chiral Ward identity ]3l|| . With 50 lattices, we only had a signal at k — 0.118, 
where we found Za = 1.210(8). This is again quite close to the pure kinetic 
result of 1.237. Dividing out the kinematic factors, this would give Za — 0.978. 

At /3 = 6.0, Ref. fH found Zy = 0.824, Z A = 1-09 for the usual clover 
action, and improved operators. 

It would be very interesting to measure the mixing of left- and right-handed 
operators (for Bk) in this action. 
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K 


^cons 


Z loc 




i (j« 

2 k 


0.118 


1.020(1) 


1.253(2) 


1.228 


1.237 


0.120 


1.032(2) 


1.185(4) 


1.148 


1.166 


0.122 


1.054(5) 


1.129(6) 


1.071 


1.098 


0.123 


1.071(6) 


1.127(8) 


1.052 


1.065 



Table 3: Vector current renormalization factors, C sw = 1.2, (3 — 5.7. 

6 Conclusions 

We have shown that via a combination of fattening the links and tuning the 
magnitude of the clover term, it is possible to optimize the chiral behavior of 
the clover lattice fermion action. The example we presented reduces the spread 
in the low energy real eigenmodes by about a factor of three in units of the 
squared pion mass, compared to the usual Wilson action. We fixed the fattening 
and varied the clover coefficient, but it is clear that a real optimization would 
involve varying both factors. 

The same procedure can be applied to the usual clover action, without fat 
links. Then, the value of the clover term C sw which minimizes the spread of 
real eigenmodes due to instantons is so large, that new kinds of configurations 
become exceptional. Their singular behavior is unrelated to fermion modes sit- 
ting on individual instantons. It is most likely related to short distance vacuum 
fluctuations. Thus this choice of action does not represent a good choice for 
improvement. 

The method of construction exposes an apparent fundamental upper limit 
for the lattice spacing in a QCD simulation of about 0.2 fm. At larger lattice 
spacing, the doubler modes and the low energy modes do not show a clean 
separation, and the mechanism of chiral symmetry breaking is qualitatively 
different than in the continuum. The Wilson and clover actions, and the one 
improved action we tested, the hypercubic action of Ref. Q , failed at lattice 
spacing 0.24 fm. This test does not apply to actions involving heavy quarks, 
for which chiral symmetry is (presumably) not important, but before any other 
improved action can be said to reproduce continuum light hadron physics at 
some lattice spacing, it should show a separation between the spectrum of its 
near low energy modes and its doubler modes. 

The "fat clover" action we have presented in this paper, as an example of 
an action with improved chiral properties, has many other nice properties as 
well. It appears to exhibit scaling of hyperfine splittings at 0.2 fm. The action 
is quite insensitive to the UV behavior of the underlying gauge field. It has 
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Z-factors for simple matrix elements which are very close to unity. And finally, 
the amount of resources required to construct propagators is halved compared 
to the usual clover action, at equivalent parameter values. It would be trivial 
to modify any existing clover code to improve it in the way we have described. 
Of course, its kinetic properties, including the dispersion relation, artifacts at 
large am q , and power law scaling violations in matrix elements, are unchanged 
from the standard clover action's. 
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